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Recently we proposed a microscopic approach to the description of the phase behaviour and critical 
phenomena in binary fluid mixtures. It was based on the method of collective variables (CV) with a 
reference system. The approach allowed us to obtain the functional of the Ginzburg-Landau- Wilson 
(GLW) Hamiltonian expressed in terms of the CV (fluctuating densities). The corresponding set of 
CV included the variable connected with the order parameter. In this paper we use this approach 
to the study of the critical behaviour of ionic fluids. For the restricted primitive model (RPM) we 
obtain the functional of the grand partition function in the phase space of the two fluctuating fields 
conjugate to the fluctuating densities. First we calculate the phase diagram of the RPM in the 
mean-field (MF) approximation and then we do this calculation taking into account the terms of 
the higher orders in the effective Hamiltonian. In the both cases the phase diagrams demonstrate 
the gas-liquid (GL) and charge ordering phase instabilities. In the latter case, the obtained value for 
c£h the GL critical temperature is in good agreement with the MC simulation data whereas the critical 

density is underestimated. The explicit expression found for the grand thermodynamic potential in 
^ . the vicinity of the GL critical point implies a classical critical behaviour of the RPM. 

?3 ' PACS numbers: 05.70.Fh, 05.70.Jk, 02.70.Rr, 64.70.Fx 
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I. INTRODUCTION 
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O , Nowdays the theory of phase transitions and critical phenomena is well developed in general. It enables one to 
predict both universal and non-universal properties for many model systems. However, a number of questions still 
remains open, among which the critical behaviour of ionic fluids is of great interest. Recent experiments have shown 
^ ■ that both Ising-like and mean-field like criticality can be observed in such systems. For reviews of the experimental 
Q\ ' and theoretical situation see Refs. Q, 0, IS El S • Numerous theoretical and computer simulation studies of the 
restricted primitive model (RPM), the simplest model for ionic systems, have not provided a clear picture of the 
thermodynamics in the critical region [1 H la ITfl ITTL 111 HI HI fl5| . 

In this paper we address the issue of the critical behaviour of the RPM using the approach proposed in 0, Il7j 
• for the binary symmetrical mixture. The theory has its origin in the approach based on a functional representation 
of a partition function in the collective variables (CV) space 0, 0] . Its particular feature is a choice of the phase 
space in which the system is considered. Among the independent variables of this space there could be the ones 
connected with the order parameters. This phase space is formed by a set of CV. Each of them is a mode of 
density fluctuations corresponding to the specificity of the model under consideration. This approach allows one to 
determine, on microscopic grounds, the explicit form of the effective GLW Hamiltonian and then to integrate the 
partition function in the neighborhood of the phase transition point taking into account the renormalization group 
symmetry. As a result, non-classical critical exponents and analytical expressions for thermodynamic functions are 



obtained [HEI] More recently this theory was developed for a binary fluid mixture [Tl ItR EH ElEl . 

The paper consists of two parts. In the first part we obtain the functional of the grand partition function (GPF) of 
the RPM given in the phase space of the two fluctuating fields: the fluctuating field hk—o conjugate to the fluctuating 
total number density /?k=o and the fluctuating field 7k conjugate to the fluctuating charge density Ck- Restricting 
our consideration to the second powers of 7k (taking into account the higher powers of /ik=o) we derive the equation 
for the chemical potential. Based on the chemical potential obtained from the linearized equation (that corresponds 
to the mean-field (MF) approximation) we calculate the spinodal curve. Its run suggests that two types of phase 
instabilities can occur in the RPM. One part of the spinodal is the gas-liquid (GL) type while another one looks 
like a A-line. We obtain the following values of the GL critical point: T* = 0.084 and r\ c = 0.005 which agree with 
the other MF theories 26]. In order to study the nature of the criticality of the RPM as well as to get the best 
estimations for its GL critical point we go beyond the above mentioned approximation taking into account the terms 
of higher orders in the effective Hamiltonian. The second part of the paper is devoted to this end. First, we consider 
the Gaussian approximation of the functional of the GPF. It yields the equation for the boundary of stability with 
respect to the charge density fluctuations. Then, applying the procedure proposed in 0,0] we obtain the expression 
for the grand thermodynamic potential in the vicinity of the GL critical point as a power series in the field /ik=o (up 
to h 4 ) conjugate to the order parameter. The expression obtained has the form of the Landau free energy. We also 
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calculate the phase diagram demonstrating both GL and charge ordering phase instabilities. The data for the GL 
critical point are T* = 0.0502 and r] c = 0.022. 



II. FUNCTIONAL REPRESENTATION OF THE GPF OF THE RPM. AN EQUATION FOR THE 
CHEMICAL POTENTIAL; THE PHASE DISGRAM IN THE MF APPROXIMATION 

The RPM consists of N = N + + iV_ hard spheres of diameter a with N + carrying charges +q and N- (= N + ) 
charges —q, in a medium of dielectric constant D. The interaction potential of the RPM has the form 

tj i \ / 00 if r < cr 

Uys(r) = | y if r > cr ' * = ±3- 

We split the potential U 1 s{r) into short- and long-range parts using the Weeks-Chandler- Andersen partition [27j . 
As a result, we have 



Vv( r ) = 




if r < a 
if r > a ' 



*7*M = 



if r < cr 



^ if r > cr 

Dr 



This simple form for $ 7< 5(r) inside the hard core changes the behaviour of the Fourier transform for large k from 
usual Colombic fc -2 to k~ 3 decay. As was shown |28(, this choice of $ 7 a(r) for r < a produces rapid convergence of 
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the series of the perturbation theory for the free energy. The Fourier transform of <fr 1 g (r) = -g- = $c(r) has the form 

(3p$ c {x) = (1) 

where j3* — jy^, (3 — f\ = § per 3 is fraction density, x = ka. 

We start with the GPF. We use the formulas obtained within the framework of the CV method for a binary 
continuous system (see ||3HJ|-flI2Jl in Appendix A). In the case of the RPM we have V(k) = 0, U(k) = 0, W(k) = 2$ c (fc) 
and = (fix,- (Mlji i s a P ar t of the chemical potential of the ith species, i — +,—). As a result, we can present the 
functional of the GPF of the RPM in the form (for details see Appendix A ) : 



3 = 3 J(dp)(dc)exp ^f3pip ~ ^ ^ c (k)c k c^ J(p,t 



(2) 



Here So is the GPF of the RS. In the case of the RPM, the RS is a one-component hard-sphere system with the 
diameter cr (potential ip~fs{f) = ip(f))- 

Pk and Ck are the CV connected with total density fluctuation modes and charge density fluctuation modes, 
respectively. p\ (pi = + /xi._)/\/2) is determined from the equation 



dlnSi _ (JV> 
9/3^1 ~ \/2* 

For the RPM J(p, c) has the same form as that for the symmetrical binary fluid 17] : 

J(p,c) = y"(d7)(dw) exp | i27r^(w k Pk + 7kCk) - «27r^ ^9Jl^ 0) w k (5k 

(-i27T) 2 1 V^/^fOl _f2) s (-l27T) 3 1 

+ 2^T*~ ^ ( 2 kCJ ~ k + * 7k7 ~ k) + 

V « k V * 



(3) 



3 
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^ (9K^ 0) cj kl cj k2 cj k3 +32n^ ) w kl 7k 2 7k 3 )^k 1 +k 2 +k3 1 



ki,k 2 ,k 3 



4! 



: — ^ (9Jl^ 0) w kl cjk 2 Wk 3 Wk4 + 69Jlf ) cJk 1 Wk 2 7k 3 7k4 + 9^1 

v2 ki,...,k 4 
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(-i2n) 5 1 ^ (o) 
X7k 1 7k 2 7k 3 7k 4 )0k 1 + ...+k4 H ^5 (9^5 ^k^k^ka 



\/2 kl ,...,l 



XWk 4 Wk 5 + 109Jt^ 2 ) Wk 1 Wk 2 Wk 3 7k4 7k5 + 5OT^ 4) Wk 1 7k 2 7k 3 7k4 7k 5 ) 

(-z2^) 6 1 ^ (o) 

X(5 kl+ ... +k5 + ^ -g 2^ W w ki u lo u k3Wk4Wk 6 Wk (1 

v2 k!,...,k 6 

+ 159Jl|j 2) cj kl cj k2 cj k 3W k4 7k 5 7k 6 + 159Jtg 4) Wk 1 ^k 2 7k3 7k4 7k 5 7k 6 ) 

X*k I+ ...+ke + •••}• ( 4 ) 

In 10} the cumulants OJln"' 1 with i„ = are connected with the nth structure factors of the RS pjjj: 

SErtW = (iV}S n . 

Structure factors S n (0) with n > 2 can be obtained from S2 (0) by means of a chain of equations for correlation 
functions (2^. Cumulants with i n ^ can be expressed in terms of SJli ' (see also formulae (4.8) in pflj): 

ot! 2 > - anw i5 an< n 4 > = 3ajii 2 2 - 2<°2 3 , 

Ml 6) = 150^3 - 30aj^ 2 4 + 16JW^ B . (5) 
First we integrate in J5J over the variables p k and w k . As a result, we have 



J2j~r h o J y (^ c ) cx p ^--^^$c(fc)c k c_ k ^ J ( c ' /l ° 



3 = aoexp | > . ^-j-I'll | /('/' ! <^p j --> /l> r (/.')r u r_ u ) ./(r. /,„). ((,) 

where 



J(c,h ) = I (dry) exp ( i2n ^ 7 k c k + ^ ^ ^ XI 7k7-k ^ OT 

J \ k Z ' V2 k n>0 



(2) K 
n+2 . 



(~^) 4 1 \^ * . *S 



4! 



X 7k x 7k 2 7k 3 7k4<5k 1+ ...+k4 X OT «+4"IT ' ( 7 ) 



_q ^ ........ . . - ^ ., , . n , 

ki,...,k 4 n>0 

^•k=o = h = ^j=- (8) 

Expressions 10-13 do not include the "field" variable /ik with 0. As one can see below, this fact will give rise to 
the Landau type free energy of the RPM in the vicinity of the GL critical point. 
Performing in ©-0 integration over Ck we obtain (up to 7 4 ) 



and condition © has the form: 



X ( 1 + V $c(fc) ^r) +L ^! 7? E 7k l7 k 2 7k 3 7k4 

ki,...,k4 



x*k 1+ ... + k4(3^-2-^ U , (10) 



where 



u n>l v ; 
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u n>2 v ' 

and formulas JSJ) are used for 2Jt^™\ 

It should be point out that this integration makes sense for a positive function In our study $c(fc) (see 0) 

can change the sign and the transition from ©-0 to (|1C)|1 has to be regarded as an algebraic formal way of deriving 
the expansion in terms of ho and 7k solely. As a result, in IjlUI) in the exponent we have the expression in terms of 
the two fluctuating fields, namely, the fluctuating field /ik=o conjugate to the fluctuating total density pk=o and the 
fluctuating field 7k (for all of k) conjugate to the fluctuating charge density Ck- 

a. Equation for the chemical potential. Restricting our consideration in to the second power of 7k and as- 
suming > (nJt^j^Tj > 0), we integrate in ©-0 over 7 k and Ck- As a result, we obtain for 



n 



V w dho 

Using and (|13[) we can write the following equation for the chemical potential 



exp(^(/io)). (13) 



where ft , .F(fto), 9 g^ o) and are given by © and JTrj-JTlJl, respectively. 

6. Phase diagram of the RPM in the MF approximation. We solve equation (|14|l in the simplest approximation. 
Neglecting terms h,Q, h$ and etc. in the right hand side of Ijl4(l and setting 

^=<<> (15) 



where 



in the left hand side of (|14fl we obtain for ^1,4. (= Mi.-) 

Mi,+ = -^a(/3), (16) 

° W = ~?i + #*cM<>- <17) 

On the other hand, we can get the same result for taking into account the terms proportional to 7 2 and 7 2 /io 

in the exponent of and setting 9Jl^ = in the final result (OJtg 2 -' is the coefficient of 7 2 /io). Therefore, we may 
say that jUi + given by l|ltj|) - (|17|) is obtained in the MF approximation. 
The full chemical potential /x+ is equal to p2| : 

where jUo,+ (= Mo,-) is the chemical potential of a one-component hard sphere system. 

The equation /3p^p" — 0, where p is the total number density, gives the spinodal of the RPM in the approximation 

considered. Using for <tc(k) and the Perkus-Yevick approximation for the RS, the equation for the spinodal curve 
can be written as 

x 2 (sina;) 2 dx tt (1 + 2?7> 2 

(18) 



(x 3 T* + 24r/sinx) 2 24r) (1 - r]) 4 

The phase diagram calculated from l|18|) is shown in Fig. 1. The obtained data for the GL critical point are 
T* ~ 0.084 (T* = 1/0*) and n c ~ 0.005. 

It is evident that in the considered approximation (corresponding to the MF theory) we get the overestimated value 
for the GL critical temperature and the underestimated value for the critical density. But, in contrast to the previous 
results, the spinodal curve changes its run (at rj ~ 0.047) and then it directs to the higher temperature. The second 
positive slop of the spinodal indicates another type of the phase instability appearing in the RPM (similarly to the 
A-line in a symmetrical binary fluid) . Below, we shall calculate the phase diagram of the RPM taking into account in 
the terms of the higher orders. 
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III. CRITIC ALITY OF THE RPM: BEYOND THE MF THEORY 

Now let us rewrite ©-(0) as 

Si = exp(J r (/i )) y"(rfc)exp ^-^-^$ c (fc)ckC-k^ 

x J {di) exp | i2vr ^ 7k c k + ^ fp X] 7k7-k9tti 0) 



k — v-s k 

2 1 /A-EYJ,^ ,_.\ , ;->-•, i 



(- ? 2tt) 2 1 ^ /a^(/i ) -(oA , (-^) 4 1 



E 7k l7 k 2 7k 3 7k 4 f 3^ - ) . kl+ ... +k4 J (19) 



X 

ki,...,k 




and the notations in (|19|l are the same as those in the previous section. 

a. Gaussian approximation. First, we restrict our attention to the second powers of variables ho and 7 k which 
corresponds to neglecting the terms proportional to /ig, 7k7— k^O) 7k7-k^0' ^o> e ^ c - After integration in 119|) over 7k 
we obtain 



= exp(^ ))n(^ 0) ) _1/2 /(d C )exp{-E^ 

k J [ k 3Jtl 



CkC-k 

(0) 



As is seen from (J^UJl, the equality 



[l + ^l>c(fc)97li 0) K. (20! 



1 + ^$ c ( fc )OT(°) = (21) 



holds at some values of the wave-vector k, temperature and density. Equation (|21|) determines the boundary of 
stability connected with the charge fluctuations (the field variable 7k is conjugate to the CV c k ): 

T* = -24^, T* = i-, (22) 
x A p* 



or 



where x* is determined from the condition 



T ;(x* ) r,) = -8r 1 ^, (23) 



tanir* = (24) 



which yields x* ~ 4.0783. Substituting x* in (|S01) we obtain the boundary of stability with respect to fluctuations of 
the local charge density 

T*(x = x*) ~ 0.28577. (25) 

A similar result (for another choice of interaction inside the hard core) was obtained in [9lll0| within the framework of 
the field-theoretical approach. The possibility of the charge-ordering transition in the continuous-space RPM model 
was discussed therein. 

It is worth noting that the RPM does not demonstrate the GL phase transition in the Gaussian approximation. 
In order to obtain the GL spinodal curve we should take into consideration the terms of the order higher than the 
second one (j 2 h, 7 2 /i 2 , etc.). 
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b. Beyond the Gaussian approximation (model if 4 )- We restrict our consideration in l|19l) to the terms of the 
fourth order. In this case 2 has the form: 

S = E exp (jF(/i )) y*(dc)exp (-^^l> c (fc)ckc_k j J {dy) exp {i2it 



k 2!V2- k 

E 7k 1 7k 2 7k 3 7k 4 (3OT(° ) -2< ) ),V ki ,... ... Ui ). ,■>(,, 

where the summation in over n is restricted to 4. 

Now we follow the programme proposed in |23| for a two-component fluid system. First, we separate the two 
types of variables: the essential variables (which include the variable connected with the order parameter) and the 
non-essential variables. Then, integrating over the non-essential variables with the Gaussian density measure, we 
construct the basic density measure (the GLW Hamiltonian) with respect to the essential variables. 

For the RPM in the vicinity of the GL critical point the variable ho (conjugate to the CV po) turns out to be the 
essential variable 01- Thus, we can present ({2l))l as 

X (l + A+±A 2 + ..) exp(-E c k c^/Tl[ 0) ) , (27) 



k 



where 



After integration in (|27|l we get 



a = ( h m { 2 ] + Km^ 



dckC- 



Si = II / = - «P E ^ , (28) 

k y'l + ^cW \n>i n -J 

where 

M n = 9JtL 0) + A93T„, 
A9Tt„ are the corrections obtained as the result of integration over variables c^: 

AWlr = A9Jl 2 = igjtf 508), 

A2Jt 3 ee 0, A9Jt 4 = 0. (29) 

and 5(/3) is given by (|T7I) . 

It is worth noting that the non-zero corrections A9Jl„ (see l|29|l ) are those which include only one sum over k. 
Next, the shift is carried out in order to eliminate the cubic term in (|28|l 



where A = —jfc- Then (gSJ has the form 



h = h + A, 



Si =Cexp I Y,Mn\ ] , (30) 
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where 

C = n . 1 :exp(<) 
22J^ U) 2 

In H30|l /To is the field conjugate to the order parameter for the GL critical point (see (1.8)-(1.10) in [2l| and (15)-(20) 
in U) 

Now we can obtain the grand thermodynamic potential of the RPM in the vicinity of the GL critical point 

n-n Q + k B T\nC = -k B T [ ^M n ^- ] . (31) 

,n>l 

The right hand side in (|31[) has the form of the Landau free energy expressed in terms of the field h conjugate to the 
order parameter. From the equation M.2 = we obtain the equation for the GL spinodal curve 



a(j3) = -2 



9Jt 2 0) OTf 



m (o) m (0)- 



c 2 smxdx _ 2 52(0) _ 5a(0) 



x 3 T* + 24r]smx 5 3 (0) 5 4 (0)' 

where S n (0) is the nth structure factor of the one-component hard-sphere system at k = 0. 

The phase diagram of the RPM is shown in Fig. 2. The curve with the maximum is the GL spinodal calculated 
using lj3~2l . The Percus-Yevick approximation is used for 52(0) (the expressions for 52(0), 53 (0) and 5 4 (0) are given in 
Appendix B). The straight line calculated by (|25|l (the Gaussian approximation) corresponds to the charge ordering 
phase transition. The GL critical point is located at T* = 0.0502 and i] c — 0.022. While the value for T* is in good 
agreement with the recent data of computer simulations 0,0] (T* ~ 0.05), the critical density is underestimated 
(?7 C ~0.04). 

We can also obtain from l|31|l (using @) the expression for 

1 Jn 2 a \P> 3 /m?<») 



,- (33) 



If we set mf ] = the expression 

(ESI reduces to lfTo|) . 

It should be pointed out that the above described scheme of integration in the vicinity of the GL critical point 
cannot be used in the region close to the line defined by (|25(l (dotted line in Fig. 2): the variables 7k appear to be 
the essential variables in this region. 



IV. CONCLUSIONS 



We use the recently developed approach in the study of the critical behaviour of ionic fluids. For the RPM we 
obtain the functional of the GPF in terms of the two fluctuating fields, namely, the field h conjugate to the total 
density fluctuations and the field 7k conjugate to the charge density fluctuations. The phase diagram of the system 
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is calculated in the MF approximation (that corresponds to setting 9Jt„ — for n > 3 in the expression for the 
chemical potential) as well as taking into account the higher powers of the field h . In the both cases the phase 
diagram demonstrates the GL and charge ordering phase instabilities. In the latter case the obtained value for the 
GL critical temperature correlates well with the MC simulation data. 

Following the scenario proposed in [Hi [T?) for a binary fluid we integrate over 7k (non-essential variables) with the 
Gaussian density measure. As a result, we obtain an explicit expression for the grand thermodynamic potential of the 
RPM in the vicinity of the GL critical point as an expansion in terms of the field conjugate to the order parameter. 
The expression implies the classical critical behaviour of the RPM. 



Appendix A 

Let us consider a classical two-component system of interacting particles consisting of N a particles of species a and 
Nb particles of species b. The system is in volume V at temperature T. 

Let us assume that an interaction in the system has a pairwise additive character. The interaction potential between 
a 7 particle at ?i and a 8 particle at may be expressed as a sum of two terms: 

where ^^(r) is a potential of a short-range repulsion that can be chosen as an interaction between two hard spheres 
<7 77 and ags- 4>is( r ) is an attractive part of the potential which dominates at large distances. 
Let us start with a grand partition function of a two-component continuous system 



£ 



£ 



N a =0 N b =0j=a,b 



n ih / (dr)exp 



7 5 ij 



(34) 



where {dT) — Yl~/dTN-,, cTV T = dr\dr\ . . .drj^ is an element of the configurational space of the 7th species; z 7 is 
the fugacity of the 7th species: 27 = exp(/3/z 7 ), p 1 = p 1 + fi^ 1 ln[(27rm 7 /? _1 ) 3 / 2 /ft. 3 ]; (3 — , &b is the Boltzmann 
constant, T is temperature; m 7 is mass of the 7th species, h is the Planck constant. /i 7 is determined from 



ains 



where (N~) is the average number of the 7th species. 

Further consideration of the problem is done in the extended phase space: in the phase space of the Cartesian 
coordinates of the particles and in the CV phase space. An interaction connected with the repulsion (potential 
tp-sir)) is considered in the space of the Cartesian coordinates of the particles. We call this two-component hard- 
spheres system a reference system (RS). The thermodynamic and structural properties of the RS are assumed to be 
known. The interaction connected with an attraction (potential </> 7d -(r) ) is considered in the CV space. The phase 
space overflow is cancelled by introduction of a Jacobian of the transition to CV. The contribution of the short-range 
forces to the long-range interaction screening is ensured by averaging this Jacobian over the RS. 

Let us introduce the grand partition function of the RS 



00 00 



Af a =0 N b =Q 7=a,6 



exp (P$N y ) 



7V 7 ! 



(dT) exp 



-?££^M ? 



7<5 ij 



(35) 



where /Xq is the chemical potential of the 7th species in the RS. 
Then the grand partition function (|34ll can be written as |22, Hj| : 



- — "0-1, 



where Ho is given in (|35|l . The part of the grand partition function which is defined in the CV phase space has the 
form of the functional integral: 



(dp)exp[(3^2fiJpo n - — ^2^24> 7 s(k)pu n p-k, S ]J(Pa,Pb)- 



(36) 



7<5 k 



9 



Here, 

1) pj is a part of the chemical potential of the 7-th species 



Mi = Mt - Mo + 57 <?Vy( fc ) 



and is determined from the equation 



dhx'Ex . , T . 



p 1 is the full chemical potential of the 7-th species; 

2)Pk,7 = Pk 7 ~ i s the collective variable of the 7-th species, the indices c and s denote the real part and the 
coefficient at the imaginary part of Pk r y, p c ~ and p£. 7 describe the value of k-th fluctuation mode of the number of 
7-th species particles. Each of p£ 7 and pt takes all the real values from —00 to +00. (dp) is a volume element of 
the CV phase space: 



( d p) = n d p°-f n d pki d p\ 



k,7 ■ 



k/0 



The prime means that the product over k is performed in the upper semispace; 
3) J(p a , Pb) is the transition Jacobian to the CV averaged on the RS [23. l23|: 



J(Pa,Pb) 



exp 



7— a 



^k,7Pk,7 



exp 



E 

n>l 



(-z27r) r 



E E M 7i-7n( k l>---) k «) I/ ki ,71 • ■ • ^k„,7„ 



71 ••O'n ki...k„ 



(37) 



where variable 1^,7 is conjugated to CV /9k,- 



,(ki,. 



. , k„ ) is the n-th cumulant connected with 



S- 



7i—7) 



(fcl 



, k n ), the n-particle partial structure factor of the RS, by means of the relation 



M~ 



,(ki, . . . ,k„) = \/N 11 . . . N ln S 11 ... ln {k\, . . . , fc„)(5 kl+ ... +lt 



where (5kH (_ kn is the Kronecker symbol. 

In general, the dependence of JVf 71 . .. 7n (ki, . . . , k„) on wave vectors ki, . . . , k„ is complicated. Hereafter we shall 
replace M 7l ... 7n (ki, . . . , k„) by their values in long-wave length limit M 7l ... 7n (0, . . . , 0). 

4) (/> 7< 5(/c) is the Fourier transform of attractive potential 4> 7 s(r). The behaviour of 4>-ys(r) in the region of the core 
r < <T 7 5 must be determined from the conditions of optimal separation of the interaction. 

We pass in ()36J) to CV /?k and Ck (according to Wk and 7k) by means of the orthogonal linear transformation: 



Pk = -7T(Pk,a + Pk,b), 



Ck 



V2 



(Pk,a — Pk,b), 



Wk = — [v-k.a + fk,b), 



As a result, we obtain for S 0, 

H = S / (dp) (dc) exp 



V2 

7k = -z-(vk,a - Vk,b)- 



PptPo + PHi Co 



W^(fc)c k c-k + U(k)pkC-k] J(p, c) 



2V 



Y}V(k)p^p^ 



(38) 



Here the following notations are introduced: 
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Pk and Ck are the CV connected with total density fluctuation modes and relative density (or concentration) 
fluctuation modes in the binary system, (dp) and (dc) are volume elements of the CV phase space: 

(dp) = dpoJJ dp^dpl, (dc) = dcoJJ efc k cfc k . 

k#0 k^O 



Functions p^ and p 1 have the form: 



-yOui.a Mi = -Fr(pi,a -Hi,b) (39) 



v(k) = (4> aa (k) + 4> bb (k) + 24> ab (k))/2 
w(k) = (4> aa (k) + 4> bb (k)-24> ab (k))/2 
u(k) = (4> aa (k)-4> bb (k))/2, 



(40) 



J(p,c)= (duj) (dj) exp i2iry^ j (vkPk + 7fc c fc) J(w,7), 



(41) 



j( W)7 ) = exp[EE t fr : E ^ } (o,...,o) 



n>l i„>0 

7ki ■ • ■7k i „^k i . 



ki...k„ 
■ w k„ 



(42) 



Index i n is used to indicate the number of variables 7k in the cumulant expansion l|42Jl . Cumulants M„ ^ are expressed 
as linear combinations of the partial cumulants M 7l ... 7ri (see (J37J) and are presented for 71, . . . , 7„ = a, b and n < 4 
in (see Appendix B in pfh 



Appendix B 



S 2 (0) 
S 3 (0) 
S 4 (0) 



(l + 2?7) 2 ' 

(1 — r/) 7 (\ — 7ij — Qf] 2 ) 

(1 + 27?) 5 ' 

(1 - ?7) 10 (1 - 30?y + 8I77 2 + 140ry 3 + 6O77 4 
(l + 2?7) 8 
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FIGURE CAPTIONS 

FIG.n The phase diagram of the RPM calculated from l|18|) (see the text for explanation). 
FIG. El The phase diagram of the RPM calculated from lE5U|l (see the text for explanation). 



13 



0.08- 



0.06- 



0.04- 



0.02- 



0.00 




i — 1 — i — 1 — i — 1 — r 

0.00 0.02 0.04 0.06 0.08 0.10 



FIG. 1: The phase diagram of the RPM calculated from 1181 (see the text for explanation). 



